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Abstract 

I We calculate the number spectrum of particles radiated during a scattering into a heat 

bath using the thermal largest-time equation and the Dyson-Schwinger equation. We 
show how one can systematically calculate ^^^^^'^ to any order using modified real time 
^ finite-temperature diagrams. Our approach is demonstrated on a simple model where 

two scalar particles scatter, within a photon-electron heat bath, into a pair of charged 
CD I particles and it is shown how to calculate the resulting changes in the number spectra of 

the photons and electrons. 

(N ■ 
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m 

^ ■ 1 Introduction 
o , 

, In recent years much theoretical effort has been invested in the understanding of rela- 
tivistic heavy ion collisions as these can create critical energy densities which are large 
enough to produce the quark-gluon plasma (the deconfined phase of quarks and gluons) 

^. m. 

^ . A natural tool for testing the quark-gluon plasma properties would be to look for the 
particle number spectrum formed when a particle decays within the plasma itself. As 
I the plasma created during heavy ion collisions is, to a very good approximation, in 
I thermodynamical equilibrium (somewhat like a "microwave oven" or a heat bath), 
one can use the whole machinery of statistical physics and QFT in order to predict the 
final plasma number spectrum. Such calculations, derived from first principles, were 
carried out by Landshoff and Taylor 1^ . 

Our aim was to find a sufficiently easy mathematical formalism allowing us to perform 
mentioned calculations to any order. Because unstable particles treated in & can not 
naturally appear in asymptotic states, we demonstrate our approach on a mathematically 
more correct (but from practical point of view less relevant) process; namely on the 
scattering of two particles inside of a heat bath. The method presented here however, 
might be applied as well to a decay itself (provided that the corresponding decay rate 
is much less than any of the characteristic energies of the process). In this paper we 
formulate the basic diagrammatic rules for the methodical perturbative calculus of plasma 
particle number spectrum '^^^J^^'' and discuss it in the simple case of a heat bath comprised 
of photons and electrons, which are for simplicity treated as scalar particles. In Section 2 
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we review the basic concepts and techniques needed from the theory of the largest-time 
equation (both for T = and T 7^ 0) and the Dyson- Schwinger equation. Rules for the 
cut diagrams at finite-temperature are derived and subsequentlv extended to the case 
when unheated fields are present. It was already pointed out in a that the thermal cut 
diagrams are virtually the Kobes-Semenoff diagrams 11 in the Keldysh formalism ^ . This 
observation will allow us to identify type 1 vertices in the real time finite-temperature 
diagrams with the uncircled vertices used in the (thermal) cut diagrams, and similarly 
type 2 vertices will be identified with the circled, cut diagram vertices. As we want to 
restrict our attention to only some particular final particle states, further restrictions on 
the possible cut diagrams must be included. We shall study these restrictions in the last 
part of Section 2. 

As we shall show in Section 3, the heat-bath particle number spectrum can be conveniently 
expressed as a fraction. Whilst it is possible to compute the denominator by means of the 
thermal cut diagrams developed in Section 2, the calculation of the numerator requires 
more care. Using the Dyson- Schwinger equation, we shall see in Section 4 that it can 
be calculated through modified thermal cut diagrams. The modification consists of the 
substitution in turn of each heat bath particle propagator by an altered one. We also 
show that there must be only one modification per diagram. From this we conclude that 
from each individual cut diagram we get n modified ones (n stands for the total number 
of heat-bath particle propagators in the diagram). Furthermore, in the case when more 
types of the heat bath particles are present, one might be only interested in the number 
spectrum of some of these. The construction of the modified cut diagrams in such cases 
follow the same procedure as in the previous situation. We find that only the propagators 
affiliated to the desired fields must be altered. 

In Section 5 the presented approach is applied to a toy model in which the quark-gluon 
plasma is simulated by sca/ar photons and electrons and we calculate the resulting changes 
in the number spectrum of the "plasma" particles. 

Finally, in Appendix A we derive, directly from the thermal Wick's theorem, the (ther- 
mal) Dyson- Schwinger equation as well as other useful functional identities valid at finite 
temperature. 

2 Basic tools 

2.1 Mean statistical value 

The central idea of thermal QFT is based on the fact that one can not take the expectation 
value of an observable A with respect to some pure state as generally all states have non- 
zero probability to be populated and consequently one must consider instead a mixture 
of states generally described by the density matrix p. The mean statistical value of A is 
then 

{A) = Tr{pA), (2.1) 
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where the trace has to be taken over a complete set of physical states. For a statistical 
system in thermodynamical equilibrium p is given by the Gibbs canonical distribution 



here Z is the partition function, H is the Hamiltonian, is the conserved charge (e.g. 
baryon or lepton number), p is the chemical potential, K = H — pN, and (3 is the inverse 
temperature: /? = 1/T (A;^ = !)• 



2.2 Largest-time equation at T=0 

An important property inherited from zero-temperature QFT is the largest-time equation 
(LTE) HI 0' il. Although the following sections will mainly hinge on the thermal LTE, it 
is instructive to start first with the zero-temperature one. The LTE at T = is a gener- 
ally valid identity which holds for any individual diagram constructed with propagators 
satisfying certain simple properties. For instance, for the scalar theory one can define the 
following rules: 
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The here means complex conjugation. In addition to these rules, each 1st type vertex 
has a factor —ig and each 2nd type vertex has a factor ig. Using this prescription, we can 
construct diagrams in configuration space. With each diagram then can be associated a 
function F{xi, . . .Xn) having all the 2nd type vertices underlined. For example, for the 
triangle diagram in Fig.|l| we have 
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F(xi,X2,X3) = {-igf iApixi - X2)iAF{xi - X3)iAF{x2 - X3) 
F{xi,X2,X3) = {-igf{ig)iA^{xi-X2)iA'^{xi-X3)iAF{x2-X3) 
F{xi,X2,X3) = {igf{-ig){-i)A*p{xi-X2)iA'^{xi-X3)iA^{x2-X3) 
F{xi,X2,X3) = {igf {-i)Ap{xi - X2) {-i)Ap{xi - X3) {-i)A*p{x2 - X3). 
etc. 

The LTE then states that for a function F{xi, . . . corresponding to some diagram 
with n vertices 

F(...,x„...) + F(...,x„...) = 0, (2.3) 

provided that XiQ is the largest time and all other underlinings in F are the same. The 
proof of Eq. (|2.3|) is based on an observation that the propagator iA p (x) can be decom- 
posed into positive and negative energy parts, i.e. 



iApix) = e{xo)iA+{x) + e{-xo)iA-{x), (2.4) 
tA^ix) = j ^^e-'^-e{±k,)5{k'-m'). (2.5) 

Incidentally, for Xjo being the largest time this directly implies 

lApi^Xj — Xi) = iA^{xj—Xi), 

iApiO) = -iA*p{0). (2.6) 

As F{. . . ,Xi, . . .) differs from F{. . . , x,, . . .) only in the propagators directly connected to 
Xi - which are equal (see Eq.( p.6| )) - and in the sign of the Xi vertex, they must mutually 
cancel. 

Summing up Eq.( p.3|) for all possible underlinings (excluding xi), we get the LTE where 
the special role of the largest time is not manifest any more, namely 

F{xi,X2,...,Xn) =0. (2.7) 

index 

The sum Y^index means summing over all possible distributions of indices 1 and 2 (or 
equivalently over all possible underlinings). The zero-temperature LTE can be easily 
reformulated for the T-matrices. Let us remind that the Feynman diagrams for the S- 
matrix (S" = 1 + iT) can be obtained by multiplying the corresponding F{xi, . . . ,x„) 
with the plane waves for the incoming and outgoing particles, and subsequently integrate 
over Xi . . .Xn- Thus, in fixed volume quantization a typical n- vertex Feynman diagram 
is given by 
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(2.8) 



Here the momenta {pj} are attached to incoming particles at the vertices {xm }, while 
momenta {qt} are attached to outgoing particles at the vertices {x„if.}. In order to 
distinguish among various functions F{xi, . . . , Xn) with the same variables Xi, . . . , Xn, we 
shall attach a subscript to each function F. For instance, the function Fi^{xi, . . . ,X4) 
corresponding to the diagram 



Pi 




1i 



contributes to {qiq2\iT\pip2) by 



4 Q-i(pi+P2)xi gi(gi+(?2)a;4 

n dx,——==——=={iAF{x, - X2))\iA{x2 - Xs))'{iA{xs - X4))^ 



similarly, the function ^24(2:1, . . . ,2:4) corresponding to the diagram 




contributes to (?i?2|«7'|piP2) by 



4 g-i(pi+P2)a;i gj(qi+?2)a;4 

TT dxi — ^= ^^ ^iApixi - X2)iAF{xi - X3){iAF{x2 - xa))^ 

X iAp{x4 — X2)iAF{x4 — X2), 



etc. 



This can be summarized as 



■ ■ ■ / n E n -1^=^ n ■■■^^n), (2.9) 

•' i=l In 3 X^^Vj^ k 



Consider now the case = |{9fc}) (let us call it |a)). From the unitarity condition: 

T — = iT^T, we get 

{a\T\a) — {a \T\a)* = t{a\T^T\a). (2.10) 
On the other hand, by construction F{xi, . . . , Xn) = F*{xi, . . . , Xn), and thus (see (^.7D) 

F(xi, ...,Xn)+F*{Xu...,Xn) = - ^ ^(^^1' ■■■^^n). (2.11) 

index' 



The prime over index in ( |2.11| ) indicates that we sum neither over diagrams with all type 



1 vertices nor diagrams with all type 2 vertices. Using ( p.9|) , and identifying with 

\{Pk}) {=\a)) we get 



{a\T\a) - {a\T\a)* = - ^ (a|T|a), (2.12) 

J , index' 

{a\T^T\a) =i ^ {a\T\a). (2.13) 



index 



Eg. ( p. 121) is the special case of the LTE for the T-matrices. The finite-temperature 



extension of ( [^.13| ) will prove crucial in Section 4. 

Owing to the 9{±ko) in A'^i^x) (see Eq.(p.5|)), energy is forced to flow only towards type 
2 vertices. From both the energy-momentum conservation in each vertex and from the 
energy flow on the external lines, a sizable class of the diagrams on the RHS's of (|2.12D- 



( p.l3| ) will be automatically zero. Particularly regions of either 1st or 2nd type vertices 
which are not connected to any external line violate the energy conservation and thus 
do not contribute (no islands of vertices), see Fig.|^. Consequently, the only surviving 
diagrams are those whose any 1st type vertex area is connected to incoming particles 
and any 2nd type vertex area is connected to outgoing ones. From historical reasons the 
border between two regions with different type of vertices is called cut and corresponding 
diagrams are called cut diagrams. 




Figure 2: An example of a cut diagram in the ip^ theory which does not contribute to the 
RHS's of ( \2.1^ )- ^.l!^ ). Arrows indicate the flow of energy. 

We have just proved a typical feature of T = QFT, namely any cut diagram is divided 
by the cut into two areas only, see Fig.|^. Eq. (|2.12|) , rewritten in terms of the cutings is 
so called cutting equation (or Cutkosky's cutting rules) H il . 
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Figure 3: Generic form of the cut diagram at the T = 0. Shadow is on the 2nd type 
vertex area. 



One point should be added. Inserting the completeness relation J2f I/) (/I = 1 into the 
LHS of (I2I3I ), we get 



J:{a\T^f){f\T\a) = ^Y^{a\T\a). (2.14) 



cuts 



It may be shown o @ that all the intermediate particles in |/) correspond to cut lines. 
This has a natural extension when {a\T^T\a) — > (alT'^VTla) with V being a projection 
operator ("P = = V^) which eliminates some of the states |/). It is easy to see that in 
such case 

{a\T^VT\a) =iJ2{a\T\a), (2.15) 

cuts 

where tilde over the J2cuts indicates that one sums over the diagrams which do not have 
the cut lines corresponding to particles removed by V. 

There is no difficulty in applying the previous results to spin-^ §0. The LTE follows 
as before: the diagram with only iSp propagators (and —ig per each vertex) plus the 
diagram with only (iSp) propagator f\ (and ig per each vertex) equals to minus the sum 
of all diagrams with one up to n — 1 the type 2 vertices {n being the total number 
of vertices). For gauge fields more care is needed. Using the Ward identities one can 
show that type 1 and type 2 vertices in (|2.12|) - (|2.13|) may be mutually connected only 



by physical particle propagators (i.e. neither through the propagators corresponding to 
particles with non-physical polarizations or Fadeev-Popov ghosts and antighosts). 

2.3 Thermal Wick's theorem (the Dyson-Schwinger equation) 

The key observation at finite temperature is that for systems of non-interacting particles 
in thermodynamical equilibrium Wick's theorem is still valid, i.e. one can decompose the 
271-point (free) thermal Green function into a product of two-point (free) thermal Green 
functions. This may be defined recursively by 

*The function iSF{x), similarly as {iAp)* (x), interchanges the role S'^ and S~. Unlike bosons, for 
fermions iSpix) is not equal to {iSF)*{x). Despite that, Eq.( 2.12| ) still holds @. 
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(r(^(xi) . . . = E ep{T{^{x,)^{xj)))^ (r( n (2.16) 



where ep is the signature of the permutation of fermion operators (= 1 for boson op- 
erators) and T is the standard time ordering symboL We shall use, from now on, the 
subscript /3, emphasizing that the thermal mean value describes a system in thermody- 
namical equilibrium (at the temperature /9~^). Note that the choice of "i" in ( |2.16| ) is 
completely arbitrary. The proof can be found for example in i^L Similarly as at 
T = 0, Wick's theorem can also be written for the (free) thermal Wightman functions 
1,0, i.e. 

{lp{xi) . . . ij{x2n))p = ^P{'^i^l)^i^j))l3 ( n ^i^k))p- (2.17) 

A particularly advantageous form of this is the so called Dyson-Schwinger equation (see 
Appendix A) which, at the T ^ 0, reads 



(G[^]^(x)F[^])^ = I dz{i^{xMz))JG[^]^^) +/ci^(^(z)^(x))^/^[^ 



5V^(z)//i"^^^'"^^^'"^^^^^\5^(z)^^'"V, 

(2.18) 

where iIj{x) is an interaction-picture field and G[. . .] and F[. . .] are functionals of ip. 
The arrc 
namely: 



The arrowed variations are defined as a formal operation satisfying two conditions. 



^ ■(^„.(x)^,(y)) = ##V^,(l/) + (-l)^V'..(x)^, (2.19) 



or 



with 



S^pniz) S^pniz) ^pniz) 

. , . ^ = 6{x-y)6mn- (2.21) 

The "p" is for bosons and 1 for fermions; subscripts m, n suggest that several types of 

fields can be generally present. Note, for bosons = ^ which we shall denote as 
For more details see Appendix A. 



2.4 Thermal largest- time equation 

The LTE (|2.13| ) can be extended to the finite-temperature case, too. Summing up in 
( ^.ltj| ) over all the eigenstates of K {= H ~ fiN) with the weight factor e~^^' {i labels 
the eigenstates), we get 
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{TT^)p = i T.{T)p. (2.22) 

index 



Let us consider the RHS of ( p.22| ) first. The corresponding thermal LTE and diagram- 
matic rules (Kobes-Semenoff rules 0) can be derived precisely the same way as at T = 
using the previous, largest-time argumentation It turns out that these rules have 

basically identical form as those in the previous Section, with an exception that now 
(0| . . . |0) {. . .)f3. Note that labeling vertices by 1 and 2 we have naturally got a dou- 
bling of the number of degrees of freedom. This is a typical feature of the real-time 
formalism in thermal QFT (here, in so called Keldysh version 0). 

We should also emphasize that it may happen some fields are not thermalized. For exam- 
ple, external particles entering a heat bath or particles describing non-physical degrees 
of freedom ¥3. Particularly, if some particles (with momenta {pj}) enter the heat bath, 
the mean statistical value of an observable A is then 



p-l3K, 

i ^ 

P{p,} = \{Pj}){{Pj}\, 
which we shall denote as {A){p.j^f^. From this easily follows the generalization of ( |2.22| ) 



{TT^){P.},^ = ^ E/^)fe},/3- (2-23) 

index 

Unlike T = 0, we find that the cut diagrams have disconnected vertex areas and no 
kinematic reasonings used in last Section can, in general, get rid of them. This is because 
the thermal part of {(p{x)ip{y))i3 Q describes the absorption of on shell particle from the 
heat bath or the emission of one into it. Thus, at T 7^ 0, there is no definite direction 
of transfer of energy from type 1 vertex to type 2 one as energy flows in both directions. 
Some cut diagrams nevertheless vanish. It is simple to see that only those diagrams 
survive in which the non-thermalized external particles "enter" a diagram via the 1st 
type vertices and "leave" it via the 2nd type ones. We might deduce this from the 
definition of {T)^p.y^p, indeed 

index index * 

Note, we get the same set of thermal cut diagrams interchanging the summation Y^^^dex' 
with J2i- It is useful to start then with Z^index' fell^K' fe})- T^i^^s is, as usual, 
described by the (T = 0) cutting rules. In the last Section we learned that the general 
structure of the corresponding cut diagrams is depicted in Fig.|^, particularly the external 
particles enter the cut diagram via type 1 vertices and leave it via type 2 ones. Multiplying 

tNote that {(p{x)<f{y))p = (: ^{xMy) :) p + {0\vixMy)\0) and (: f{xMy) ■)p ^ J ^ fB{ko)6{e - 
^2y-~tk(x-v) ^ with /B(fco) = (e'^l'^ol - 1)-^ 
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each diagram (with the external particles in the state |i; {pj})) with the prefactor ^ ^ 
and summing subsequently over i, we again retrieve the thermal cut diagrams, though 
now it becomes evident that the particles {pj} enter such diagram only via type 1 vertices 
and move off only through type 2 ones, since the summation of the (T = 0) cut diagrams 
from which it was derived does not touch lines corresponding to unheated particles. Note, 
the latter analysis naturally explains why the unheated particles obey the (T = 0) LTE 
diagrammatic rules even in the thermal diagrams 

Another vanishing comes from kinematic reasons. Namely three-leg vertices with all 
the on shell particles (1-2 lines) can not conserve energy- momentum and consequently 
the whole cut diagram is zero. As an illustration let us consider all the non-vanishing, 
topologically equivalent cut diagrams of given type involved in a three-loop contribution 
to*E..wmp.,/3 (see Fig.|). Q 



— 













2^ 


~\2 








y2^s 











a) b) c) d) 







y^c — 




yr^C — 




>(V — ' 






/7\ 






/7 iV 


y2\ 







e) f) g) h) 

Figure 4: An example of non-vanishing cut diagrams at the T 7^ 0. The heat-bath consists 
of two different particles. External particles are not thermalized. 



Let us stress one more point. In contrast with T = 0, at finite temperature the cut itself 
neither is unique nor defines topologically equivalent areas, see Fig.^, only the number of 
crossed legs is, by definition, invariant. 







Figure 5: The cut diagram from Fig)^ c) demonstrates that the cut can be defined in many 
ways but the number of crossed lines is still the same. 



This ambiguity shows that the concept of the cut is not very useful at finite temperature 
and in the following we shall refrain from using it. 

In Chapter 4 it will prove useful to have an analogy of ( |2.23| ) for {T^VT)^. Here V is the 
projection operator defined as 

■'■Let us emphasize that originaUy we had the 64 possible cut diagrams. 
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V = Y.\a;j){a;j\, (2.25) 

j 

where "j" denotes the physical states for the heat-bath particles and "a" labels the 
physical states for the outgoing, non-thermalized particles. Let us deal with (T^PT)^. 
Using ( |2.15| ), we acquire 

{T^VT), = zY^^ Y.{l\m- (2.26) 

' index' 

Interchanging the summations, we finally arrive at 

{T^VT)^ = tJ2{T)p, (2.27) 

index 

where tilde over the J2index' means that we are restricted to consider the cut diagrams, 
with only (l-2)-particle lines corresponding to the "a" and "/'particles (i.e {0\ip{x)(p{y)\0) 
and {ip{x)ip{y)) /3, respectively). The extension of Eq.( p.27|) to the case where some ex- 
ternal, non-thermalized particles {pk} are present is obvious, and reads 

{T^'PT){,,},p = ^t.{T)^,,},p. (2.28) 

index 

Finally, let us note that using the LTE, one may extend the previous treatment to var- 
ious Green functions. The LTE for Green's functions is then a useful starting point for 
dispersion relations, see e.g. fllUl. 



3 Heat-Bath particle number spectrum: 
general framework 

The cutting equation (|2.28|) can be fruitfully used for both the partition function Z 
and the heat-bath particle number spectrum '^^^^"•'^ calculations. To see that, let us 
for simplicity assume that two particles (say $i,$2) scatter inside a heat bath. We 
are interested in the heat-bath number spectrum after two different particles (say 0i, ^2) 
appear in the final state. Except for the condition that the external particles are different 
from the heat bath ones, no additional assumption about their nature is needed at this 
stage. 

The initial density matrix pi (i.e. the density matrix describing the physical situation 
before we introduce the particles $i(pi), $2(^2) into the oven) can be written as 

A = Zr'Y.^~'''''\r,PuP2){r,Pi,P2\, (3.1) 
j 

where "j" denotes the set of occupation numbers for the heat bath particles. A long time 
after the scattering the final density matrix pf reads 
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pf 



(3.2) 



here V is the projection operator projecting out all the non-heat-bath final states except 
of 4>i{qi) , (p2{(l2) ones. The S'-matrix in (|3.2|) is defined in a standard way: 5=1 + iT. 



The Zf in ( |3.2| ) must be different from Zj as otherwise pf would not be normalized to 
unity. In order that pf satisfy the normalization condition Tr{pf) = 1, one finds 



Zf = T.^-^'''^r:Pi.P2\S^VS\3-p,,P2) = {S^VS),,,,,p Z, = {T^VT),,,,^p Z,. (3.3) 

The key point is that we have used in ( p. 3D the T - matrix because the initial state 
|$i(pi), $2(^2)) is, by definition, different from the final one |0((?i), 02('?2)) and conse- 
quently VS can be replaced by iVT . This allows us to calculate Z^ using directly the 
diagrammatic technique outlined in the preceding Section. 

From ( p.l|) and (|3.2| ) one can directly read off that the number spectrum of the heat bath 
particles is: 



dw J (27r)3 ^ 



and consequently 



_ 1^2 _ ^2^ {T^'Pa}i^;^)aii^;^)T)p,p,,f3 



where we have used the completeness relation for the final states |/) and [P; a^a] = 0. 
The subscript "/" denotes which type of heat-bath particles we are interested in. In the 
following the index will be mostly suppressed. 



4 Modified cut diagrams 



To proceed further with (|3.4|) and ( |3.5| ), we expand the T- matrix in terms of time-ordered 



interaction-picture fields, i.e. 

TVP] = ^ i dxi. . . j rfx„a„(xi . . . Xn)T{'4){xi) . . . ipixn)). (4.1) 

Here ip represents a heat-bath field in the interaction picture. Other fields (i.e. 0,0 and 
$) are included in the a„ |. An extension of ( [4.1|) to the case where different heat- 
bath fields are present is natural. Employing ( |4.1D in {T'^VT)p-^p^^p, one can readily see 

^ When Fermi fields are involved, we have, for the sake of compactness, included in the argument of 
■0 the space-time coordinate, the Dirac index, and a discrete index which distinguishes ipa from -0^. 
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that this factorizes out in each term of the expansion a pure thermal mean value (. . .)/?. 
The general structure of each such thermal mean value is: {Gm['ip]Fn[ip]) p, where . .] 
and Gm[- ■ ■] are the operators with "n" chronological and "m" anti-chronological time 
ordered (heat-bath) fields, respectively. Analogous factorization is true in the expansion 
of {T'^V^ a'^ aVT) pj^p2^i3. The only difference is that the pure thermal mean value has 
the form {Gm['ip]a'^(iFn[ip]) 13 instead []. In case when various heat-bath fields are present, 
m = mi + 7712 + . . . + mn, with "m;" denoting the number of the heat-bath fields of l-th 
type. 

Applying the Dyson-Schwinger equation to {Gm[4']0'^ CLFn[ip]) /3 twice and summing over 
"n" and "m", we get cheaply the following expression (c.f. also ( |A.ll )) 



{T^VajaiT) 



PlP2,/3 ~ 



A similar decomposition for {T'^VT) p^p^^^ would not be very useful (cf. ( [A.18| )); instead 
we define {{T^VT)')p ^P2^I3 having the same expansion as {T'^VT) p^p^^p except for the 
an{- ■ .)Valj^{. . .) are replaced by «„(. . .)Val^{. . .) "-'+™' . In this formahsm {{T^VT) )pip2,/3 
is 



dxdy I T'^S 



2 



, (4.3) 



with the T being the anti-chronological ordering symbol. Comparing ( [4.3| ) with ( [A.19| ), 



we can interpret the RHS of ([4.3|) as a sum over all possible distributions of one line 
(corresponding to tpi) inside of each given (T 7^ !) cut diagram constructed out of 

^Remember that V — V ®V" — g2)('J'i, 92! ®Yl,]\i){3\- Here V' — behaves as an 

identity in the subspace of heat-bath states. 
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{T'^VT)p^p^^p. As ( f4.3| ) has precisely the same diagrammatical structure as {T''Va'^aT)p-^p^^i3- 
{a^ a) ^{T'^VT) p^p^^p (cf.( [4.2|) ), it shows that in order to compute the numerator of '^^^^'■'^-'^ = 

'^^^duj^'^^ — '^^^j^'^''' f\ one can simply modify the usual {T^'PT)p-^p^^p cut diagrams by the 
following one-line replacements (cf. (|3.4|) ). 



(i) For neutral scalar bosons: 

+ {ip{x)a{k;uj))f3{a\k;uj)ip{y))fs} 

S{e-m'){fB{u;){fB{u;) + l) 



(27r)3 

+ S+ih -co){l + fBico)) - 5-{k, + u)fB{uj)]e~'^^^-y\ (4.4) 

where fsi^) is the Bose-Einstein distribution: /b(i^) = ^i3M_i - Term 0{—kQ)fB{uj) de- 
scribes the absorption of a heat-bath particle, so reduces the number spectrum, that is 
why the negative sign appears in front of it. Analogously, 

(T(v,(xMy))), ^ j ^^6\uj'-k'-m'){{a\k-uj)^{x)),{a{^-uj)^{y)), 

+ {a{k]uj)^{x))p{a\k]uj)^{y))p] 

d^k 



^^^^-5{k' - m2)(l + fs{uj))fB{u)e'''^-y^ 
X {5+{ko-uj) + 5~{k^ + uo)). (4.5) 



Similarly, for A{N) one needs the following replacements (cf. ( p.5D ) 

{^{xMy))p - /(05(^'-^'){/b(^.)(/b(c..) + i) 

+ e{ko){i + fBiook)) - e{-ko)fB{ook)} e-'^^^-y\ 

{TicfixMy))), l0^,5{e-m')il + fBM)fBMe-^'^^-y\ 



(4.6) 



with uJk = V^^-m?. 



llHere^M^ = /^<5+(c.2-k2-m2)(at(a;,k)a(c.,k))^, (cf. Q). 
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(ii) For Dirac fermions: 

The Dirac field is comprised of two different types of excitations (mutually connected via 
charge conjugation), so the corresponding number operator N{u!) — Ni,{u) + Nci{u) with 



a=l,2'' l^^J 



a=l,2' 



Thus, the one-line replacements needed for '^^(^bi^^)) ^^^^ 



a=l,2 



(27 

- {M^)b^i^;uj))f,{bUk;u;)i^M)p} 



X {(1 - /^(a;)) - Mu;){l - ^(a;))}e-^'=M, 

(4.7) 

where /^(a;) is the Fermi-Dirac distribution: fF{(-o) = and 



-— ^ 5+(u;2 - - m') {(6„(k;u;)V'p(x));3(&i(k;u;)V.(y))/3 
-(6t(k;a;)^,(a;));,(6«(k;a;)^,(|/))4 



x/p(u;)(l-/^(a;))e-'=M. 



(4.8) 



Correspondingly, for A(iVft) we need 
— f d^k 

X {(1 - /^(o;)) - /^(a;)(l - ^(a;)}e-'=(--^) 

(4.9) 
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For the d-tjpe excitations the prescription is very similar, actually, in order to get 
the following substitutions must be performed in (|4^ )- (|4.9| ): 9{ko) 9(—ko), 
/f (1 - fp) and /i -fi. 



(iii) For gauge fields in the axial temporal gauge = 0): 

The temporal gauge is generally incorporated in the gauge fixing sector of the Lagrangian 
and particularly 

Cfix = -^iAo)'';a-^0. (4.10) 

The principal advantage of the axial gauges arises from the decoupling the F-P ghosts 
in the theory. This statement is of course trivial in QED as any linear gauge (both 
for covariant and noncovariant case) brings this decoupling automatically OJ. Particular 
advantage of the temporal gauge comes from an elimination of non-physical scalar photons 
from the very beginning. 

Let us decompose a gauge field Ai, z = 1, 2, 3 into the transverse and longitudial part, i.e. 
Ai = Af + Af with 

AJ = (5., - ^) A, and Af = ^A„ (4.11) 
and use the sum over gauge-particle polarizations 

j:s['\k)sf\k) = 5,,-^, (4.12) 

A=l 

with e^^\k) being polarization vectors, then for "we get the following one-line 

replacements 

{Aj{x)Aj{y))^ ^E^y^^^(^'-k^-"^')(^n^)4(k;a^)),(aA(k;.;)4(y)), 
+ {Aj{x)a,{k;u))^{ai{k;u)Aj{y))^} 

did 



)(Eq.(4.4)) 



r H^k 

{nAj{x)Aj{y)))^ - ^E^y (^^^(^'-k^-"^'){(^f(^)4(k;^)),(An2/)aA(k;^)), 

+ {Aj{x)ax{k;uj))p{Aj{x)al{k;uj))p 

= (6.,-^]iEq.(^)). (4.13) 



Q2 
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The replacements needed for A{N ) can be concisely expressed as 



{■■■), 



5i 



^) (Eq.(4.6)) 



As for the longitudial (non-physical) degrees of freedom, it is obvious that 



(4.14) 



(4.15) 



Eqs.( |4.4| )-( |4.14| ) can be most easily derived in the finite volume limit, e.g. for a scalar 
field we reformulate ^p{x) as 



At 



iErt — lWrl^ 



(4.16) 



rescaling the annihilation and creation operators by defining a{k) = ^lE^y Aj, in such a 
way that [Ak, A^ = 6,,' (so that (^1^,.) = d.^'feih)), while / t^t ^ ^ Ek- 



The replacements ([4.4| )- (^4.14|) are meant in the following sense: firstly one constructs 
all the T 7^ diagrams for {T'^VT)p^p^^f^, using the LTE ( p.28[ ) and the rules mentioned 
therein. In order to calculate the numerator of (|3.4| ) or (|3.5| ) we simply replace (using 
corresponding prescriptions) one heat bath particle line in each cut diagram and this 
replacement must sum for all the possible heat-bath particle lines in the diagram. If more 
types of heat-bath particles are present, we replace only those lines which correspond to 
particles whose number spectrum we want to compute (see Fig.|^). 











— ~4r^ 










j7 iV 




y2 

























/2V 







Figure 6: The numerator of o-nd ^3.d[) can he calculated using the modified cut 

diagrams for {T'^VT)p^p^^p. As an example we depict all the possible contributions to 
the numerator derived from the cut diagram on Fig}^ c). The wavy lines and thin lines 
describe the heat-bath particles. The crossed lines denote the substituted propagators, in 
this case we wish to calculate the thin-line particle number spectrum. 



The terms in the replacements ( |4.4| )-( ^rT4D have a direct physical interpretation. The 
f{uJk) and (1 -|- (— 1)^/(0;^)) can be viewied as the absorption and emission of the heat- 
bath particles respectively O. The term /(u;fc)(l -|- (—1)^/(0;^)) describes the fluctuations 
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of the heat bath particles. This is because for the non-interacting heat-bath particles 
{{nk — {uk) is)'^) 13 = /(co'fc)(l + {—iyf{iJk)). The substituted propagators can be therefore 
schematically depicted as 



X y ~ fluctuations 

X \y y ~ fluctuations + emissions + absorptions 

2 ^> 1 

X y ~ fluctuations 

Collecting all the contributions from emissions, absorptions and fluctuations separately, 
one can schematically write 

(i(A^(cu)) J d(^N(^U!))i ^ pemission _|_ pahsorption _|_ J?/'""^^;^^ ^4 

where, for instance for neutral scalar bosons 



r^emission 



PlP2,/3 

Using ( [4. 51) , it is easy to write down the analogous expressions for the p"-'^^°^p^'^°'^ and 
pfiuc^ To the lowest perturbative order, the form ( [4.17 ) was obtained by Landshoff and 
Taylor @. 

5 Model process 
5.1 Basic assumptions 

To illustrate the modified cut diagram technique, we shall restrict ourselves to a toy 
model, namely to a scattering of two neutral scalar particles $ within a photon-electron 
heat bath, with a pair of scalar charged particles 0, left as a final product. Both initial 
and final particles are supposed to be unheated. We further assume that the heat bath 
photons A and electrons \1/ are scalars, i.e. the heat-bath Hamiltonian has form 

2 

W = lid^Af-^J^A^ (5.1) 

On the other hand ifj„ entering in the T - matrix reads ifj„ = + (eA + Yyl^)^\E'"l' + 

{eA + ^A^)(f)(f)^ . It is usually argued Ull that the interacting pieces in H^^ can be 



dropped provided that tj —oo and t/ ^ oo. In the following we accept this omission 
as it allows us to use safely Wick's theorem (|2.16|) and Dyson-Schwinger equation ( p.l8|) 
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5.2 Calculations 

We can now compute an order-e^ contribution to the The evaluation of the 

'^^^^^dw^^^ straightforward. In Fig.|^ we list all the modified cut diagrams contributing 
to an order-e^. 









% 












/ 1 \ 






































■~-x-' 






Figure 7: The modified cut diagrams involved in an order-e^ contribution to the photon 
number spectrum. Dashed lines: photons. Solid lines: (p, 0^ particles. Bold lines: $ 
particles. 



Note that diagrams b) and c) are topologically identical. Similarly, diagrams e), f), h), i) 
and j) should be taken with combinatorial factor 2 (corresponding diagrams for insertion 
of a heat-bath particle line in the bottom, solid line are omited). Of course, diagram 
g) vanishes for kinematic reasons. Let us emphasize that it is necessary to give sense 
to graphs e), h), i) and j) as these suffer with the pinch singularity; the muon-particle 
propagator {p\.2 — m?)~^ has to be evaluated at its pole because of the presence of an on- 
shell line (1-2 line) with the same momenta. Some regularization is obviously necessary. 
Using the formal identity 



(5.2) 



we discover that the unwanted mutually cancel between e) and h) diagrams (similarly 
for i) and j) diagrams). An alternative (but lengthier)) way of dealing with the latter 
pinch singularity; i.e. switching off the interaction with a heat bath in the remote past 
and future, is discussed in @. Using the prescriptions ( |4.4| )-(|475|) and evaluating the 
integrals (note, we should attach to each digram the factor ^ coming from a Taylor 
expansion of the T-matrix), we are left with (c.f. Eq.( [4.17D ): 
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and 



;^emiss^on _j_ jpabsorption 

tAV 

{TVT%^p,,f,V8ujp^ujp,{27r)^ 



d^k 6{k^ - mi) 6{ko - u) 



X Jd%d% S^qf-ml) S^ql-ml) 



X 



{K,{l + f{Lu))6\-Q + qi + q2 + k) 



-K2fBH 5\-Q + qi + q2-k)\ (5.3) 



(TPTt)p^P2,/3^8^Pi^P2(2vr)5 i 

X jd%d% SHqf-ml) 6+{ql-ml) 

X {5\~Q + qi + q2 + k)K^ + 5\-Q + gi + - k)K2 
- 28\-Q + qi + q2)K^] 

tAV/.M(l + /.M) r _ _ 

< J d\d% 6\-Q + qi + q2) 



+ 



2gi/i; — mi 2qik + y 2 ^ o 



+ ^ 52)} (5.4) 



2 / \ 2 



with /^i = ^2gifc+m2 + 2g2fc+m2 j ' -^2 = (^2gifc-m2 + l^^F^ j ; -^3 = (2gifc-m2)(2g2fc+m2) 

and Q = + P2- We have dropped the prescription in the propagators since adding/ 
subtracting an on-shell momenta gi;2 to/from an on-shell momenta k we can not fulfil 
the condition {k ± gi;2)^ = "o^^- As it is usual, we have assumed that our interaction is 
enclosed in a 'time' and volume box {t and V respectively). The relevant (i.e. order-e°) 
term for {TVT'^)p^p2,p reads 



{TVT^)p,p,,p = gy^ !i[{27ry I "^'^^"^'^^ " ""^^ ^""^^^ ~ ""^^ '^'^"^ + + 



64yu;p,cc;p2|(5|vr 



g2-4m2 (5.5) 



Eqs. (|5.3| ) and (|5.4| ) are analogous to the result obtained in for the decay. 
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One can perform the similar calculations for the electron number spectrum — • 
The task is now, however, tougher. The major difference in comparison with the photon 
number spectrum is that the lowest order in e (keeping A^) is e^. This brings richer 
diagrammatic structure then in the photon case. In Fig.^ we list all the generating 
thermal diagrams contributing to an order-e^. 




d) e) f) 



Figure 8: The generating thermal diagrams involved in an order-e^ contribution to the 
photon number spectrum. Dashed lines: photons. Thin lines: (p, (f)^ particles. Bold lines: 
$ particles. Half-bold lines: electrons. 

It is easy to see that out of these 6 generating thermal diagrams we get 41 non- vanishing 
cut diagrams (10 from a); 9 from b); 4 from c); 8 from d); 4 from e) and 6 from f).) The 
actual electron number spectrum calculations are thus rather long and not extremely 
rewarding for our unphysical toy model, so we refrain from doing them. 



6 Conclusions 

In this paper we have formulated a systematic method for studying the heat-bath particle 
number spectrum using modified cut diagrams. In particular, for the quark-gluon plasma 
in thermodynamical equilibrium our approach should be useful as an effective alternative 
to the Landshoff and Taylor ^ approach. The method used in S (i.e. to start from 
first principles) suffers from the lack of a systematic computational approach for higher 
orders in coupling constants. One of the corner stones of our formalism is the largest-time 
equation (LTE). We have shown how the zero-temperature LTE can be extended to finite 
temperature. During the course of this analysis, we have emphasized some important 
properties of the finite-temperature extension which are worth mentioning. Firstly, many 
of kinematic rules valid for zero-temperature diagrams can not be directly used in the 
finite-temperature ones. This is because the emission or absorption of heat-bath particles 
make it impossible to fix some particular direction to a diagrammatic line. It turns out 
that one finds more diagrams then one used to have at zero temperature. The most 
important reductions of the diagrams have been proved. The rather complicated struc- 
ture of the finite-temperature diagrams brings into play another complication: uncutable 
diagrams. It is well known that at zero temperature one can always make only one cut in 
each cut diagram (this can be viewed as a consequence of the unitarity condition). This 
is not true however at finite temperature. We have found it as useful to start fully with 
the LTE analysis which is in terms of type 1 and type 2 vertices. This language allows us 
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to construct systematically all the cut diagrams. We have refrained from an explicit use 
of the cuts in finite-temperature diagrams as those are ambiguous and therefore rather 
obscure the analysis. 

The second, rather technical, corner stone is the thermal Dyson- Schwinger equation. We 
have developed a formalism of the arrowed variations acting directly on field operators. 
This provides an elegant technique for dealing in a practical fashion with expectation 
values (both thermal and vacuum) whenever functions or functionals of fields admit the 
decomposition ([A.l[). The merit of the Dyson- Schwinger equation is that it allows us to 



rewrite an expectation value of some functional of field in terms of expectation values 
of less complicated functionals. Some illustrations of this and further thermal functional 
identities are derived in Appendix A. 

When we have studied the heat-bath particle number spectrum, we applied the Dyson- 
Schwinger equation both to numerator and denominator of corresponding expression. 
The results were almost the same. The simple modification of one propagator rendered 
both equal. We could reflect this on a diagrammatical level very easily as the denominator 
was fully expressible in terms of thermal cut diagrams. Our final rule for the heat-bath 
particle spectrum is 



dA{Niu)) {T^VT) 



M 



with T being the T-matrix, V being the projection operator onto final states, pi,p2 being 
the momenta of particles in the initial state, f3 being the inverse temperature and M being 
abbreviation for the modified digrams. Modification of the cut diagrams consist of the 
substitution in turn of each heat-bath particle line by an altered one. This substitution 
must be done in each cut diagram. Replacement must be only one per modified diagram. 
Our approach is demonstrated on a simple model where two scalar particles ("pions") 
scatter, within a photon-electron heat bath, into a pair of charged particles ( "muon" and 
"antimuon" ) and it is shown how to calculate the resulting changes in the number spectra 
of the photons and electrons. 
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A Appendix 

A.l Functional formalism: general background 

Eq.( |2.18D gives us an alternative definition of the Dyson- Schwinger equation in terms of 
the "functional derivation" ^J^^^ . Let us first show that ( |2.18|) is consistent with Wick's 
theorem ( p.l6|) -( p.l7 ). To be specific, let us consider an ensemble of non-interacting 
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particles in thermodynamical equilibrium. In order to keep the work transparent, we 
shall suppress all the internal indices. There is no difficulty whatsoever in reintroducing 
the necessary details. Let us first realize that for any (well behaved) functional the 
following Taylor's expansion holds 
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= Y1 / dxi... dxna^{xi . . .Xn)i'ixi) . . .^(Xn), (A.l) 

The same is true if if) is an operator instead. In the latter case the a"'{. . .) are not 
generally symmetric in the x's When Fermi fields are involved, we might, for the sake 
of compactness, include in the argument of ip the space-time coordinate, the Dirac index, 
and a discrete index which distinguishes ipa from tp^. In the latter case J dx ^ J2 I dx, 
whera summation runs over the discrete indices. With this convention, the expansion 
( |A.1| ) holds even for the Fermi fields. An extension of ( |A.1|) to the case where different 



fields are present is natural. Particularly important is the case when is a field in the 
interaction picture, using Wick's theorem and decomposition ( |A.1| ) one can then write 



{G[i:]ij{x)F[i:])p = 

= E (/ dxj (/ dyY a"(xi . . . Xn)nyi ■■■ym)(^{f[ H^k)^ H^) Wvt 

/ I- \m rn I m 

+ Hi dy) /3™(yi...y„)E(±iy"'w^M^/))HG[^]n^(?/. 



k')] . (A.2) 



with (/ dxY = J dxi . . . J dxn- The ' — ' stands for fermions and '+' for bosons. On the 
other hand, using the formal prescriptions (|2.19|) and ( p.21|) for one can read 



dz{ij{x)ij{z))JG[^] 



6F\ 



Y.[ dy) r{yi---ym) dz{^{x)ij{z))pY.(±iy-^5{z-yi)I^G[ij]\{^{yk>)) 

m J I \ k'^l I p 

Y.( dy) r{yi---yrn)Y.i±ir'{i^{xMyi))AG\i;]X{i^{yk')) • (A.3) 



Similar expression holds for \ dzl(^{x)^{z)) p (^f^ifi^M^ • Putting latter two together 

we get precisely (|A.2|) . This confirms the validity of (|2.18|) . It is easy to persuade oneself 
that exactly the same sort of arguments leads to 



**If X = the a" may also contain derivations working on the various fields. In this paper 

we rule out such a case from our reasonings. 
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{T{ij{x)Fmp = I dz{T{^{x)m)), { T I 1 ) (A.5) 



6^{z) 



5^Ij{z) 



dzmz)^{x)), [j^'^im))^ , (A.6) 



etc. 



with T being either the chronological or anti-chronological time ordering symbol. At 
this stage it is important to realize that from the definition of j^^^ directly follows that 

( ' — ' holds for bosons and '+' for fermions). Indeed, 



\-StP(x) ' <5i/)(i/)J=F 



6tlj{x)6ilj{y) 



n=2 i<j 



(A.7) 

Similarly [^^; = 0. Analogously we might prove 



6ip{x)6ilj{y) 6ip{x)6il){y) 
and 



5i){x)5i){y) 5i){x)5i){y) ^ ' 5i){x) 5^{y) 

6^{y) 6^{x) ^^^Silj{x)6^{yy ^ ' ' 

The "p" is for bosons and 1 for fermions. With (|2.18|) and (A. 4) - (A.6) one can easily 
construct more complicated expectation values. For example, using ( p.l8[ ) and (A. 4) we 
get 
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(A.IO) 

■5 



Similarly, using ( |2.18|) and (anti-)commutativity of the arrowed 777-^, we get 



{G[^]^{x)^{y)F[^])p 
dzidz2 



2 mxMz,))f,{ij{y)^{z2))p + {-irmxMz2))pmymzi))p) 

+ J dz,dz2mz,)i^{x))p{i^{y)i^{z2))p + {-lY{i^{x)i^{z2))p{i^{z,)i^{y))p) 
\5^{zr)5ij{z2)l p 

+ {ij{x)^{y))p{G[ij]F[^])p. (A.ll) 

We could proceed further having still higher powers of fields and variations. However, 
there is a quite interesting generalization in case when we have (anti-)time ordered oper- 
ators. Let us have F[iIj\ = T{F[ip]), in this case 

n ^ i=l 

n=l ^ ' ij mj^ij 

= /</..</..(T(V.(.0^W..))>.(^-|3^L) +(f[0]>,. (A,12) 



where F[ip] differs from Flip] in the replacement . .) — * ° ^"^ {n starts from 1 !). In 
comparison with (A.4)-(A.ll), the a^{. . .) (i.e. the pure T = contribution) does matter 
here. Note that . .) generally involves non-heat bath fields with corresponding space- 
time integrations. Similar extension is true if F[%p\ = Tc{F[iIj]), where 7^ is the time path 
ordering symbol. In that case 
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dzidz2{Tc{ij{zi)tp{z2))) p 



/3 

6m 



Ic \5'ilj{z2)5ip{zi) / ^ 

+ {m)p. (A.13) 

with ^(^dx = J^dt Jydx and = 5c{x — y) Wick's theorem for the 7c-oriented 
product of fields has an obvious form 

{Tc{i^{x,)...i^{x2n)))^ = Y.^p{Tc{i^{x^)i>{x,)))p {TciXH^ix,))) ^. (A.14) 

This can be directly derived from Wick's theorem ( p.l7| ), realizing that 

Tc{ip{xi) . . . ipixm)) = Y^epOcitp,, . . . , tpj-^{xp,) . . . ipixpj, (A.15) 

p 

where P refers to the permutation of the indices and 6c{ti, . . . , tm) being a contour step 



function Iz2f defined as 



^ ^ , f 1 (ti, . . . , t„ are 7^-oriented along C) ,^ 

- ■ I (otherwise) ^ ' ^ 



Particularly important is the Keldysh-Schwinger path o see Fig.l 
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tj - i e 



Figure 9: The Keldysh-Schwinger time path. 

In the latter case 



{F[^])/, = I dz^dZ2{r{ij{zMz2)))f,' ^^^^ 



Ci \6'ip{z2)6ip{zi) 



/3 



^^A contour (5-function Sc{x — y) is defined as J„dzSc{z — z )f{z) = f{z ), see HI 



c 
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+ / dzidz2{T{'4>{zi)^{z2))) p ' 



C2 \6'ip{z2)6ip{zi 



+ 2 dzi dz2{i){zi)ip{z2))f3 



C2 Jci \6iIj{z2)6iIj{zi) 



/3 



+ im)^- (A.17) 

Application to the product G['^]F['iIj] with = Tc,iF[ip]) and = Tc^{G[ip]) is 
straightforward and reads 



{G[i:]F[ij])p = I dz,dz2{T{i;{z,)i;{z2)))J ^f"^^^' 



5'iij{z2)5i){ziy 



dz^dZ2(J{,lj{z,)ij{z2)))p ( G[^]- ^'^^"^^ 



6ip{z2)6^lj{zi) 



G[iP]6 6F[^] 



+ 2 dzidz2{ij{zi)ij{z2))p . I 
J \ di){zi) di{j{z2) I p 

+ (G[0]F[0])^, (A.18) 

where the overhning indicates that we work with " '^"n^rJ" '^ instead of a^{. . . .)? 

we have also abbreviated J^^ dz to / dz /^^ dz. We should also emphasize that used 
in ( |A.18|) is S{x — y) rather than 6c{x — y). 

In Eq.( |4.3| ) it has been used the inverted version of ( [A.18| ), namely 



Here has the coefficients a" (. . .)/?™(. • .)^^ instead of . • •)• Note, 

that the . . .) does not contributes and thus we do not have any pure T = 
contributions. Eq. ( |A.19| ) has a natural interpretation. Whilst the LHS tells us, that from 
each thermal diagram (constructed out of {G[^]F[%1}\) j^) with ^^-^ internal heat-bath 
particle lines we must take n + m identical copies, the RHS says, that this is virtually 
because we sum over all possible distributions of one heat-bath particle line inside of the 
given diagram. The pictorial expression of (|A.19|) is depicted in Fig.|l^ . 
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X (number of lines) = ^ 

Ik 




Figure 10: Diagrammatic equivalent of Eq.( \A.l{^ . The cut separates areas constructed 
out of Flip] and Glip]. 
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